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Visceral  leishmaniasis  is  a  potentially  fatal  dis¬ 
ease  alTecting  large  populations  in  many  pans  of 
the  world.  In  Latin  America  alone,  some  1 .6  mil¬ 
lion  people  are  considered  at  risk,  and  almost 
200,000  cases  occur  annually.'  The  causative 
agent  is  Leishmania  donovani  chayasi.  which  is 
transmitted  by  the  Phlebotomine  sand  fly  Lui- 
zomyia  longipalpis.  Dense  populations  of  Lit. 
iongipalpis  arc  frequently  associated  with  human 
habitation  in  many  of  the  drier  regions  of  Latin 
America.  The  species  has  been  found  in  all  known 
foci  of  visceral  leishmaniasis  and  is  a  proven 
vector  in  some  of  them.-  '  More  recently,  L.  d. 
chagasi  has  also  been  isolated  from  persons  with 
a  nonulcerative  form  of  cutaneous  leishmaniasis 
in  Central  America."  ’  Lutzomyia  longipalpis  was 
the  predominant  anthropophilic  sand  fly  in  these 
foci  as  well. 

The  geographic  distribution  of  Lu.  longipalpis 
extends  from  Mexico  to  southern  Brazil  and  con¬ 
sequently  covers  a  variety  of  habitats.  This  large 
range  is  subdivided  by  mountain  ranges  and  zones 
of  unsuitable  climate  that  apparently  pose  for¬ 
midable  barriers  to  migration.  Phlebotimine  sand 
flies  are  generally  poor  fliers  with  movement  re¬ 
stricted  to  short,  flight-assisted  hopping.  Flight 
ranges  for  several  neotropical  species  in  the  genus 
Lutzomyia  have  been  estimated  not  to  exceed 
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100  meters  in  a  24-hr  period.'  Wind-assisted  dis¬ 
persal  has  been  shown  not  to  be  significant. '  Geo¬ 
graphic  isolation  among  populations  of  a  sp>ecies 
with  low  vagility  can  promote  genetic  divergence 
resulting  from  genetic  drift  and/or  selection  for 
adaptation  to  regional  habitats.  Eventually  post- 
zygotic  or  prezygotic  reproductive  isolating 
mechanisms  (RIM)  may  e\  oh  e.  If  reproductive 
isolation  occurs,  such  populations  become  in¬ 
dependent  gene  pools,  each  representing  a  spe¬ 
cies  that  follows  a  distinct  evolutionary  trajec¬ 
tory.  Ward  and  others  have  suggested  that  Lu. 
longipalpis  may  be  a  species  complex."'  Their 
conclusions  were  based  on  a  series  of  attempted 
matings  between  populations  that  differed  in  the 
number  of  pale  spots  on  the  abdomens  of  males 
(one  pair  of  spots  versus  two  pairs).  They  did 
not  find  sperm  in  female  spermathecae  and  sug¬ 
gested  the  existence  of  pre-RlM.  In  the  few  cases 
in  which  hybrids  were  produced,  both  sexes  were 
sexually  fenile.  .Although  their  results  are  inter¬ 
esting,  the  failure  of  mating  in  caged  populations 
cannot  be  taken_as  proof  of  pre-RIM  in  nature. 
A  study  of  the  genetic  structure  of  a  single  pop¬ 
ulation  of  Lu.  longipalpis  from  Bolivia,  consist¬ 
ing  of  a  mixture  of  individuals  that  varied  in 
body  size,  revealed  no  evidence  of  significant 
genetic  divergence  among  different  body  size 
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groups."  The  purpose  of  the  work  we  report  here 
was  to  estimate  the  level  of  genetic  divergence 
among  laboratory  populations  of  Lu.  lonf^ipalpis 
that  originated  from  three  widely  separated  lo¬ 
calities  and  to  determine  if  posi-RlM  exist  that 
may  be  involved  in  the  maintenance  ofdiversitv 
in  nature. 

MATERIALS  AND  METHODS 

Sand  fly  colonies 

The  Brazil  strain  of  Lu.  longipalpis  originated 
in  Lapinha  caves  near  Belo  Horizonte.  Minas 
Gerais  (approximate  coordinates  20°0'S.  44”0'E). 
The  Colombia  colony  was  derived  from  flics  col¬ 
lected  by  R  B.  Tesh  near  Melgar,  Tolima  de¬ 


partment  (4“1  1  N.  ■’4'’!  8'\V)  in  1989.  The  Costa 
Rica  colony  was  started  with  Hies  collected  by  us 
near  Libena.  Guanacaste  province  (10‘’37'N. 
85°26’W)  in  1991.  The  locations  of  collection 
sues  arc  shown  in  Figure  1 .  .All  colonies  were 
maintained  as  described  by  Modi  and  Tesh.  • 
The  Colombia  and  Costa  Rica  colonics  were 
started  with  the  progeny  of  more  than  100  field 
collected  female  Hies. 

Isozyme  analysis 

The  whole  bodies  of  individual  adult  flies  were 
homogenized  in  7  ^1  of  distilled  water.  Homog¬ 
enates  were  applied  to  horizontal  12.5%  (w  v) 
starch  gels  for  electrophoresis  using  siai  Jard 
procedures. ' '  Five  different  buffer  svstems  were 
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Figure  1.  Localities  of  the  populations  that  pro¬ 
vided  sand  flies  for  colonies, 

used  to  maximize  electrophoretic  separation  of 
enzymes.  These  are  given,  along  with  references, 
in  Table  1.  Enzyme-speciftc  histochemical  stain¬ 
ing  procedures  were  used  to  visualize  bands  on 
gels.'*  Discrete  zones  of  staining  activity  on 
gels  were  assumed  to  be  controlled  by  single  loci 
coding  for  specific  enzyme  products.  The  prod¬ 
ucts  of  27  putative  enzyme  loci  provided  ade¬ 
quate  resolution  to  study  variability.  Loci  coding 
for  the  same  enzyme  were  designated  numeri¬ 
cally  in  sequence  from  most  cathodal  to  most 
anodal.  Alleles  were  scored  on  the  basis  of  the 
distance  bands  migrated  through  the  gel  with  ref¬ 
erence  to  the  most  common  allele,  which  was 
given  the  r/  value  of  1 .00. 

Hybridization  c.\peruncnts 

Virgin  females  (6  hr  or  less  post  pupal  emer¬ 
gence)  of  one  colony  were  placed  in  a  cage  with 
males  from  another  colon} .  Two  days  later,  flics 
were  fed  blood  and  four  days  thereafter  were 
transferred  to  an  oviposition  container.  .All  in- 
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tcrcolony  crosses  and  reciprocals  were  achieved. 
Crosses  between  individuals  from  the  same  col¬ 
ony  were  likewise  conducted  as  controls.  Eggs 
were  counted,  the  larval  hatch  was  estimated, 
and  the  mean  number  of  F,  adults  per  female 
was  calculated.  Sterility  in  hybrid  males  was  de¬ 
termined  by  microscopic  examination  of  testes. 
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The  tesics  and  distal  portion  of  the  vasa  defer- 
emia  were  dissected  and  transferred  to  a  small 
drop  of  saline  on  a  microscope  slide  and  a  cover 
slip  was  added.  The  gross  appearance  of  the  testes 
was  noted.  Gentle  pressure  was  applied  to  allow 
examination  of  the  contents  of  testes.  The  prep¬ 
aration  was  examined  under  difi'erential  inter¬ 
ference  contrast  illumination.  Sperm  viability  was 
determined  based  on  total  numbers  and  unifor¬ 
mity  of  morphology  as  compared  with  parental 
males. 


HESI  LTS 

Isozyme  survey 

Eleven  of  the  27  loci  (40.7%)  studied  were 
monomorphic  throughout  the  sample.  Three  of 
the  1 6  polymorphic  loci  were  diallelic.  and  man- 
nose-6-phosphate  isomerase  was  the  most  poly¬ 
morphic  locus,  being  coded  for  by  a  total  of  five 
alleles.  The  observed  mean  r  SEM  hetero/\- 
gosity  per  individual  (H)  throughout  the  sample 


CO  BZ 


CR 


B 


[1  234  56789  101112131415 


lW  3  4  5  6  7  8  9"  10 1112 1 


CO  CR 


BZ 


BZ  CO  CR 

bS  D 

Figure  2.  Gels  illustrating  the  phenotspes  of  certain  diagnostic  loci  for  species  in  the  Luizomyiu  i<>uyipalpis 
complex.  A.  malic  enzyme  (ME).  Lanes  1-.^  and  7-9  arc  ME  " '  homo/ygoies  and  lanes  4-ti  are  ME '  homo/\goies. 

B,  mannose  phosphate  isomerase  (MPl).  Lancs  1-5  are  MPI'”'  homozygoies.  lanes  6-8  are  MPL^'  homo/ygotes. 
lanes  9  and  10  are  MPI'-’"/MPl'-‘'  hctcrozygoies.  lane  11  is  an  MPI' "  homo/ygote.  lane  12  is  an  MPL'  ho¬ 
mozygote.  lane  13  is  an  MPP"  homozygote,  lane  14  is  an  MPI'"'  homo/ygote.  and  lane  1 5  is  an  MPP  homo/ygote. 

C,  glutamate  oxaloacetate  transaminase- 1  (GOT-1).  Lanes  1-10  are  GOTI homozygotes  and  lanes  1 1-1 5  are 
GOT-L'*  homozygotes.  D,  peptidase- 1  (PEP-1).  Lanes  1-4  are  PEP-l"”*  homo/ygotes.  lanes  5-9  are  PEP-L  “ 
homozygotes,  lanes  10-12  arc  PEP-l'"'  homo/ygotes.  and  lane  13  is  a  PEP- 1 PEP-1 "  '  hetero/ygote  (the 
products  of  the  pepiidase-2  locus  can  be  seen  as  faint  bands  towards  the  top  of  the  gel).  CO  C  olombia;  BZ 
Brazil.  CR  =  Costa  Rica. 
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Table  4 

Nei  unbiased  genetic  distances  tabow  diagonal)  and 
genetic  similarities  theh'is  diagonah.  amc  v  three 
I  utzomyia  longipalpis  laboratory  colonies 


<Jolon> 

Bra/F 

Colombia 

c  osia  Kicu 

Brazil 

- 

0.333  - 

0.255 

Colombia 

0.717 

- 

0.121 

Costa  Rica 

0.775 

0.886 

- 

was  0.057  ±  0.018  (Table  2).  Each  of  the  !''• 
polymorphic  loci  carried  alleles  that  v  e  unique 
to  a  single  colony.  Seven  loci  possessed  unique 
alleles  which  were  fixed,  or  nearly  sc  in  c>ne  col¬ 
ony  (Tabic  3  and  Figure  2).  Consequently,  the 
levels  of  genetic  divergence  a  'ng  the  three  -ol- 
onies.  as  measured  by  genetic  distance  (D)'“,  were 
very  high  (Table  4).  The  diagnostic  values  of  loci 
were  calculated  using  the  method  of  Ayala  and 
Powell;-^’  a  locus  was  considered  diagnostic  if  it 
could  be  used  to  correctly  assign  a  fly  to  the 
population  from  which  it  came  with  a  probability 
greater  than  99%.  Using  this  criterion,  there  were 
six  loci  that  were  diagnostic  for  the  Brazil  and 
Colombia  colonies,  five  loci  diagnostic  for  the 
Brazil  and  Costa  Rica  colonies,  and  two  loci  di¬ 
agnostic  for  the  Colombia  and  Costa  Rica  col¬ 
onies  (Table  5).  Using  multilocus  genotypes,  cor¬ 
rect  identifications  can  be  achieved  with  a 
probability  exceeding  99.99%  for  any  combi¬ 
nation. 

Hybridization  experiments 

Table  6  summarizes  the  results  of  hybridiza¬ 
tion  experiments.  Among  the  control  crosses,  each 
female  laid  an  average  of  2  i  .65  eggs  (not  counted 
for  the  Co'ombia  colony).  The  mean  numbc'"  of 
eggs  laid  among  females  outorossed  to  males  of 
different  strains  was  29.76.  Survival  of  F,  gen¬ 
erations  to  tiie  adult  stage  was  lowest  in  the  Brazil 
colony  (22.6%)  and  was  44.0%  in  the  Costa  Rica 
colony.  The  mean  percent  survival  to  the  adult 


stage  among  the  progeny  of  live  intercolony 
crosses  was  57.7%.  which  was  substantially  high¬ 
er  than  that  observed  in  the  control  crosses.  In 
outcrossing  colony  strains,  one  might  expect  in¬ 
creased  vigor  in  the  F,  resulting  from  hetcrosi,. 
The  reproductive  systems  of  progeny  males  ifom 
control  crosses  appeared  normal  and  cont,imed 
nur  'ious  mature  spermatozoa  (Figure  3A).  The 
gross  appearance  of  testes  and  vasa  deferentia  m 
intc'colony  hybrids  did  not  appear  different  from 
ti  ,jse  of  the  progeny  of  control  matings.  How¬ 
ever,  their  testes  contained  only  few-  normal  sperm 
a  1  .iumerous  abnormal  ones  (Figure  3B).  Only 
in  the  cross  bet  we-’u  Colombia  females  and  C  osta 
Rica  males  were  some  hybrids  observed  that  had 
testes  containing  norma!  spermaiuzi  i.  A  sample 
of  F,  males  from  this  cross  were  eleciropho’"‘sed 
and  the  gels  were  stained  for  enzyme  systems 
the*  were  diagnostic  for  the  Colombia  and  Costa 
Rica  colonies.  The  results  confirmed  that  a  per- 
,age  of  these  males  w-ere  in  fact  derived  from 
Colombia  females  that  were  not  virgin  when  used 
in  the  cross,  'f  nose  males  that  had  normal  sperm 
were  probably  the  result  of  this  contaminatiop 

DISCUSSION 

We  initially  analyzed  oui  data  by  pooling  aii 
individual  genotypes  and  treating  them  as  a  sin¬ 
gle  population.  The  level  of  genetic  polymor¬ 
phism  in  the  overall  sample  studied  was  high. 
Fiity-nine  percent  of  the  27  loci  s’ud'ed  carried 
two  or  more  alleles,  however,  the  mean  hetero¬ 
zygosity  was  relativc'v  low  (5.7%.  Table  2).  The 
reason  for  thi  ispanty  was  tnat  alleles  were  not 
distributed  randomly  among  the  three  colon^ 
populations.  The  degree  of  genetic  structuring 
among  Cvilony  populations  was  profound:  c\  cry 
one  of  the  lo  polymorphic  loci  carried  at  least 
one  allele  that  occurred  in  only  a  single  colony. 
W’hen  the  data  were  partitioned  to  determine 
genetic  variability  forecch  individual  colony .  the 
mean  heterozygosity  averaged  for  the  ilnce  coi- 


Table  5 

Diagnostic  'alues  of  seven  loci  for  distinguishing  colonies  or  Lutzomyia  longipalp-'s* 


Companson  ALO  OOT-l  Mi  l  MFM  PEP  !  POM  S  rOH 

Brazil  versus  Colombia  l.OOO  0.999  1.000  0.999  1.000  0.945t  1.000 

Brazil  versus  Costa  Rica  t  0.999  1.000  0.999  t  1000  1.000 

Colombia  versus  Costa  Rica  1.000  t  t  t  t  1.000  0.--86+ 

•  ALD  —  aldolase.  GOT •  I  ~  giutamaie  oxaloacclate  iransammasc-  i;  MF-l  =  malic  en/sme-i .  MPI  =•  ma  , nose-6  phosphate  jsomerase,  PtP- 1 
peptidase*  I :  PCM  *  D^''^sphoglucomutase,  SODH  ~  sorbitol  dehydrogenase 
t  Not  dtagnostic. 
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Table  6 

Results  of  laboratory  hybridization  bei  veen  pairs  of  three  laboratory  colonies  rt/'Luizorriyia  longipaipis* 


Female  ^  male 

Total  no. 
hloodfcd 

No  of  females 
oviposited 

Mean  no  of 
eggs  female 

Mean  no  of  F 
adiiMs  Jci  ale 

F  sperni 

S  33hi]li>  ^ 

Control 

BZ  X  3Z 

48 

35 

26.5 

6.0 

- 

CO  X  CO 

NC 

.NC' 

NC 

NC 

CR  X  CR 

20 

17 

16.8 

7  4 

- 

Experimental 

BZ  X  CO 

24 

(6 

45.1 

2''.g 

CO  X  BZ 

23 

16 

39.8 

24,0 

BZ  X  CR 

26 

20 

33.7 

24.0 

CR  X  BZ 

48 

4! 

15,8 

6.5 

CO  X  CR 

NC 

20 

NC 

10.5 

- 

CR  X  CO 

58 

50 

14.4 

7.8 

- 

*  BZ  *■  Brazil:  CO  Colombia.  NC  =  not  cour  ‘cd.  C  K  -  C  osia  Kica 
t  *  =  viable;  -  =  inviablo 


onies  was  6,2%.  esseniially  the  same  as  for  the 
overall  sample.  However,  the  proportion  of  poly- 
morphic  loci  decreased  from  0.593  in  the  overall 
sample  to  0.296  on  average  in  the  three  colonies. 

Since  the  level  of  genetic  di'/ergence  between 
colonies  can  be  best  summarized  by  calculating 
pairwise  values  for  genetic  distance,  we  used  the 
procedure  of  Nei (Table  4).  Estimates  of  genetic 
divergence  were  quite  high.  The  levels  of  genetic 
distance  we  observed  were  in  fact  equivalent  to 
those  reported  among  closely  related  species  in 
the  genus  Luizomyia  (these  included  compari¬ 
sons  between  fu.  to^ensendi.  Lit.  spindrassa.  Lit. 
youngi.  Lit.  longiflocosa.  Lit.  qita.siiow  nscndi.  and 
Lu.  sauroida).-'  In  the  same  study  Kreutzcr  and 
others  reported  that  the  genetic  distance  aver¬ 
aged  for  comparisons  between  Lit.  lonyiilocoso. 
Lu.  quajiiown.sen.di.  and  Lii.  sauroida  was  1)  ^ 
0.010.  leading  them  to  suggest  that  these  tasa 
may  be  conspecific.  The  values  we  observed  were 
an  order  of  magnitude  higher  than  this,  suggest¬ 
ing  that  the  level  of  genetic  divergence  between 
pairs  of  our  three  colonies  more  closely  resemble 
the  relationships  between  species  in  this  genus, 
rather  than  local  populations  within  a  single  spe¬ 
cies. 

The  comparisons  we  describe  here  are  among 
laboratory  populations,  and  it  is  well  known  that 
laboratory  colonies  are  affected  by  forces  that 
may  have  a  significant  impact  on  their  genetic- 
makeup.  How  much  of  the  genetic  divergence  we 
observed  may  be  an  artifact  of  the  colonization 
process?  The  most  common  change  in  colony 
genetic  makeup  is  decreased  heterogeneity.’’  This 
may  -esult  from  genetic  drift  and/or  selection.  If 
we  consider  two  colonies  of  the  same  species. 


loss  of  heterogeneity  will  significantlv  impact  ge¬ 
netic  distance  omy  if  an  allele  at  a  polymorphic- 
locus  approaches  fixation  n  one  colony  and  an 
alternate  allele  approaches  fixation  in  the  other 
In  our  study,  this  weld  have  had  to  have  hap¬ 
pened  at  eight  different  loci  (Table  3).  This  is 
highly  unlikely,  considering  that  two  of  the  col¬ 
onies  were  newly  established  using  eggs  pooled 
from  more  than  100  wild-caught  females.  The 
only  information  available  on  genetic  heteroge¬ 
neity  in  natural  populations  of  Lu.  iongipalpis 
comes  from  a  single  population  in  Bolivia,  where 
observed  heterozygosity  was  measured  at  0.03"^. 
this  value  was  lower  than  that  for  any  of  our 
colony  populations  (Tabic  2).  Thus,  if  this  value 
IS  typical,  heterozygosity  was  not  likely  to  haw 
decreased  significantly  during  colonization. 

.Anoiher  phenomenon  that  may  aflcci  the  ge¬ 
netic  makeup  ol  colony  populations  is  the  exi'-- 
tcnce  of  balanced  poly  morphisms,  which  result 
in  excessively  high  heterozygosity  at  those  loci 
affected.  '  E.xcess  vc  heterozygosity  at  a  few  ioci 
may  obscure  otherwise  lower  than  normal  values 
when  averaged  over  all  loci.  Balanced  polymor- 
phismsare  the  result  ol  linkage  to  recessive  lethal 
genes,  in  which  ease  only  hcierozygotes  survive. 

.An  analysis  of  goodness  of  fit  to  Hardy-Wcin- 
berg  genotypic  frequencies  was  performed  for  the 
three  colonies.  No  fignificant  excess  of  hetero- 
zygotes  was  observed.  We  believe  that  the  high 
level  of  divergence  between  colonies  cannot  be 
adequately  explained  by  independent  coloniza¬ 
tion.  We  therefore  conclude  that  this  divergence 
reflects  the  situation  in  nature,  and  that  our  re¬ 
sults  provide  potential  genetic  markers  that 
should  be  useful  for  identification  in  the  field. 
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Figure  3.  A.  appearance  of  normal,  mauire  sperm  o'i  Liazninvia  ionaipalpis.  B,  the  lypicai  appearaiurc  o! 
the  abnormal  sperm  produced  by  imcrcolony  f  hybrid  males,  in  this  case  the  product  of  a  cross  between  a 
Brazil  female  and  a  C  olombia  male.  (Magnification  -  ■  2.000.) 


even  in  situations  where  the  three  species  occur 
in  sympatrv.  It  follows  that  some  mechanism 
must  exist  that  permits  high  levels  of  genetic 
divergence  to  persist  among  natural  populations. 
Given  the  extent  of  the  observed  genetic  diver¬ 
gence.  it  is  likely  that  this  mechanism  severely 
limits  the  amount  of  gene  flow  among  difi'erert 
populations.  In  addition  to  geographic  isolation. 


post-RlM  would  produce  such  an  cfl'cci  and  arc 
amenable  to  analysis  in  laboratory  situations. 

In  hybridization  experiments,  which  included 
all  pairwise  matings  and  reciprocal  crosses,  we 
observed  sterility  in  male  progeny  from  all  in¬ 
tercolony  crosses  (  Figure  313).  In  each  of  the  con¬ 
trol  (intracolony)  crosses,  males  were  sexually 
fertile  (Figure  3A).  The  status  of  hybrid  females 
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was  not  assessed.  In  only  one  case  did  some  hy¬ 
brid  males  appear  to  be  fenile,  but  this  was  prob¬ 
ably  due  to  experimental  error,  resulting  from 
failure  to  separate  the  sexes  from  the  stock  colony 
before  they  began  mating  among  themselves.  This 
was  evidenced  by  the  presence  of  homozygotes 
for  maternal  allozymes  in  some  F,  progeny.  It  is 
not  likely  that  this  mating  incompatibility  are.  e 
in  colony  because  the  Colombia  colony  was  es¬ 
tablished  less  than  two  years  prior  to  the  time 
these  experiments  were  conducted  and  the  Costa 
Rica  colon\  was  in  the  second  to  sixth  genera¬ 
tion.  Hybrid  sterility  is  a  common  sequela  of 
interspecific  hybridization  and  is  generally  con¬ 
sidered  as  an  indicator  of  species  status  for  the 
two  parental  forms  producing  sterile  progeny.'' 
Our  results  indicate  that  Lu.  longipalpis  is  net  a 
single  species,  but  exists  in  nature  as  a  complex 
of  at  least  three  sibling  species. 

Genetic  di%  ergence  among  vector  species  can 
include  genes  that  affect  vectorial  capacity.  The 
best  known  examples  come  from  closely  related 
anopheline  mosquitoes  that  \  ar>’  dramatically  in 
their  imponance  as  vectors  of  malaria.  While  the 
overwhelming  majority  of  L.  z/.  chagasi  infec¬ 
tions  in  Latin  America  result  in  life-threatening 
visceral  leishmaniasis,  our  newly  characterized 
Lu.  longipa/pis  species  from  Costa  Rica  transmit 
genetically  identical  L.  d.  chagasi  that  cause  cu¬ 
taneous  lesions  and  do  not  visccralizc.''  ‘  Prelim- 
inar\  results  indicate  that  tnc  saliva  ol  this  spe¬ 
cies  docs  noi  contain  a  sasodilatory  peptide 
(Warburg  and  others,  unpublished  data),  that 
was  prcviousU  identified  m  the  saliva  ol  the  Bra¬ 
zilian  species.  ■’  Because  saiivar\  giaiu.!  issaies 
from  Lu.  longipalpts  (Brazil)  have  been  sliowii 
to  enhance  the  development  of  cutaneous  leish¬ 
maniasis  infections,-’  polymorphisms  in  salivary 
composition  may  be  an  important  factor  deter¬ 
mining  the  epidemiology  and  perhajis  natholog\ 
of  the  leishmaniases. 
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